Drude weight, plasmon dispersion, and pseudospin response in doped graphene sheets 
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Plasmons in ordinary electron liquids are collective excitations whose long- wavelength limit is rigid 
center-of-mass motion with a dispersion relation that is, as a consequence of Galileian invariance, 
unrenormalized by many-body effects. The long- wavelength plasmon frequency is related by the 
f-sum rule to the integral of the conductivity over the electron-liquid's Drude peak, implying that 
transport properties also tend not to have important electron-electron interaction renormalizations. 
In this article we demonstrate that the plasmon frequency and Drude weight of the electron liquid in 
a doped graphene sheet, which is described by a massless Dirac Hamiltonian and not invariant under 
ordinary Galileian boosts, are strongly renormalized even in the long- wavelength limit. This effect 
is not captured by the Random Phase Approximation (RPA), commonly used to describe electron 
fluids. It is due primarily to non-local inter-band exchange interactions, which, as we show, reduce 
both the plasmon frequency and the Drude weight relative to the RPA value. Our predictions can 
be checked using inelastic light scattering or infrared spectroscopy. 



INTRODUCTION 

The first theory of classical collective electron density 
oscillations in ionized gases by Tonks and Langmuir pQ 
in the 1920's helped initiate the field of plasma physics. 
The theory of collective electron density oscillations in 
metals, quantum in this case because of higher electron 
densities, was developed by Bohm and Pines [2j [3] in 
the 1950's and stands as a similarly pioneering contribu- 
tion to many-electron physics. Bohm and Pines coined 
the term plasmon to describe quantized density oscil- 
lations. Today plasmonics is a very active subfield of 
optoelectronics [4] [5] , whose aim is to exploit plasmon 
properties in order to compress infrared electromagnetic 
waves to the nanometer scale of modern electronic de- 
vices. This wide importance of plasmons across different 
fields of basic and applied physics follows from the ubiq- 
uity of charged particles and from the strength of the 
long-range Coulomb interaction. 

The physical origin of plasmons is very simple. When 
electrons in a plasma move to screen a charge inhomo- 
geneity, they tend to overshoot the mark. They are 
then pulled back toward the charge disturbance and over- 
shoot again, setting up a weakly damped oscillation. The 
restoring force responsible for the oscillation is the av- 
erage self-consistent field created by all the electrons. 
Because of the long-range nature of the Coulomb in- 
teraction, the frequency of oscillations ou p \(q) tends to 
be high and is given in the long wavelength limit by 
Upi(q — > 0) = nq 2 V q /m where n is the electron density, 
m is the bare electron mass in vacuum, and V q is the 
Fourier transform of the Coulomb interaction. This sim- 
ple explicit plasmon energy expression is exact because 
long- wavelength plasmons involve rigid motion of the en- 
tire plasma which is independent of the complex exchange 
and correlation effects that dress [6 the motion of an indi- 
vidual electron. The exact plasmon frequency expression 



is correctly captured by the celebrated RPA [2j [3] [6] , but 
also by rigorous arguments [7 in which the selection of 
a particular center-of-mass position breaks the system's 
Galilean invariance and plasmon excitations play the role 
of Goldstone bosons. In two-dimensional (2D) systems 
V q = 27re 2 /q so that co p \(q — > 0) = ^2iine 2 q/m. 

Electrons in a solid, unlike electrons in a plasma or 
electrons with a jellium model [6] background, experience 
a periodic external potential created by the ions which 
breaks translational invariance and hence also Galilean 
invariance. Solid state effects can lead in general to a 
renormalization of the plasmon frequency, or even to the 
absence of sharp plasmonic excitations. In semiconduc- 
tors and semimetals, however, electron waves can be de- 
scribed at super- atomic length scales using kp theory [8], 
which is based on an expansion of the crystal's Bloch 
Hamiltonian around band extrema. In the simplest case, 
for example for the conduction band of common cubic 
semiconductors, this device leads us back to a Galileian- 
invariant parabolic band continuum model with isolated 
electron energy E c (p) = p 2 /(2m\^). The crystal back- 
ground for electron waves appears only via the replace- 
ment of the bare electron mass by an effective band mass 
mt>. It is this type of k • p Galilean invariant interact- 
ing electron model, valid for many semiconductor and 
semiconductor hetero junction systems, which has been 
of greatest interest in solids. The absence of electron- 
electron interaction corrections to plasmon frequencies 
at very long wavelengths in these systems has been am- 
ply demonstrated experimentally by means of inelastic 
light scattering [9j [10] . 

The situation turns out to be quite different in 
graphene - a monolayer of carbon atoms tightly packed 
in a 2D honeycomb lattice [TTJ [12j [13] , which has engen- 
dered a great deal of interest because of the new physics 
it exhibits and because of its potential as a new mate- 
rial for electronic technology. The agent responsible for 
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many of the interesting electronic properties of graphene 
sheets is the bipartite nature of its honeycomb lattice. 
The two inequivalent sites in the unit cell of this lattice 
are analogous to the two spin orientations of a spin- 1/2 
particle along the +z and — z directions (the z axis being 
perpendicular to the graphene plane). This observation 
opens the way to an elegant description of electrons in 
graphene as particles endowed with a pseudospin degree- 
of-freedom [TTJ [12j [13] (in addition to the regular spin 
degree-of- freedom which plays a passive role here). When 
k • p theory is applied to graphene it leads to a new type 
of electron fluid model, one with separate Dirac-Weyl 
Hamiltonians for electron waves centered in momentum 
space on one of two honeycomb lattice Brillouin-zone cor- 
ners: Hb = va p. Here v is the bare electron velocity, p 
is the k • p momentum, and cr is the pseudospin operator 
constructed with two Pauli matrices {<r% i = x,y}, which 
act on the sublattice pseudospin degree-of-freedom. It 
follows that the energy eigenstates for a given p have 
pseudospins oriented either parallel (upper band) or an- 
tiparallel (lower band) to p. Physically, the orientation 
of the pseudospin determines the relative amplitude and 
the relative phase of electron waves on the two distinct 
graphene sublattices. 

The feature of graphene that is ultimately responsible 
for the large many-body effects on the plasmon disper- 
sion and the Drude weight is broken Galilean invariance. 
What happens is that the oriented pseudospins provide 
an "ether" against which a global boost of the momenta 
becomes detectable. This is explained in detail in the 
caption of Fig. [I] 

From Fig. [T we can also see why the plasmon fre- 
quency in graphene is so strongly affected by exchange 
and correlation. In a plasmon mode the region of oc- 
cupied states (Fermi circle) oscillates back and forth in 
momentum space under the action of the self-induced 
electrostatic field. In graphene however, this oscillatory 
motion is inevitably coupled with an oscillatory motion 
of the pseudospins. Since exchange interactions depend 
on the relative orientation of pseudospins they contribute 
to plasmon kinetic energy and renormalize the plasmon 
frequency even at leading order in q. 

In what follows we present a many-body theory of this 
subtle pseudospin coupling effect and discuss the main 
implications of our findings for theories of charge trans- 
port and collective excitations in doped graphene sheets. 



GRAPHENE DIRAC MODEL 

Graphene's honeycomb lattice has two-atoms per unit 
cell and its tt- valence band and 7r* -conduction band touch 
at two inequivalent points, K and K' \ in the honey- 
comb lattice Brillouin-zone. The energy bands near e.g. 
the K point are described at low energies by the spin- 
independent massless Dirac Hamiltonian Hb introduced 
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FIG. 1: Breakdown of Galileian invariance in 
graphene. Panel la) shows the occupied electronic states 
in the upper band of graphene in the ground state. Notice 
that every state is characterized by a value of momentum (the 
origin of the arrow) and a pseudospin orientation (the direc- 
tion of the arrow). Panel lb) shows the occupied states after 
a Galilean boost. An observer riding along with the boost 
would clearly see that the orientation of the pseudospins, rel- 
ative to the center of the occupied region has changed. It 
looks like the pseudospins are subjected to a "pseudomagnetic 
field" that causes them to tilt towards the +x direction. The 
appearance of this pseudomagnetic field is the signature of 
broken Galilean invariance. In contrast, in a Galilean invari- 
ant system [Panels lc) and Id)] the energy eigenstates are 
characterized by momentum only: an observer riding along 
with the boost would not see any change in the character of 
the occupied states. 



above. Electron-electron interactions in graphene are de- 
scribed by the usual non-relativistic Coulomb Hamilto- 
nian He , which is controlled by the 2D Fourier transform 
of the Coulomb interaction, V q = 2ite 2 /(eq) with e an ef- 
fective average dielectric constant. 

Electron carriers with density n can be induced in 
graphene by purely electrostatic means, creating a cir- 
cular 2D Fermi surface in the conduction band with 
a Fermi radius ftp, which is proportional to ^Jn[lA\. 
The model described by H = 7Yd + requires an ul- 
traviolet wavevector cut-off, /c max , which should be as- 
signed a value corresponding to the wavevector range 
over which TLb describes graphene's tt bands. This cor- 
responds to taking /c max ~ l/ a o where ao ~ 1.42 A is 
the carbon-carbon distance. This model is useful when 
& max is much larger than fcp- I n this low-energy descrip- 
tion, the many-body properties of doped graphene sheets 
depend p~5j [16] on the dimensionless fine- structure cou- 
pling constant a ee = e 2 / (ehv) (which is defined as the 
ratio between the Coulomb energy scale e 2 k^/e and the 
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kinetic energy scale hvkp) and on density via the ultra- 
violet cut-off A = /c max //cF. The fine-structure constant 
a ee can be tuned experimentally by changing the dielec- 
tric environment surrounding the graphene flake [T7J [18] . 
The ultraviolet cut-off A varies from ~ 10 for a very high- 
density graphene system with n ~ 10 14 cm -2 to ~ 100 for 
a density n ~ 10 12 cm -2 just large enough to screen out 
unintended [19] inhomogeneities. From now on Planck's 
constant h divided by 2tt will be set equal to unity, h = 1. 

The collective (plasmon) modes of the system can be 
found by solving the following equation [6] , 

1-^^(9,^ = (1) 

where XppiQ, °°) is the so-called proper [20] density-density 
response function. In the q — » limit of interest here we 
can neglect the distinction between the proper and the 
full causal response function XppiQi 00 )- We show below 
that 

v 2 q 2 

lim lim 5fte x pp (q, u) = A — r (2) 

where A is a density-dependent constant which has the 
value Ao = gsF/(4:7rv 2 ) for a ee — > 0. Here g = g s g v = 4 
accounts for spin and valley degeneracy and Sp = vkp is 
the Fermi energy. Note the order of limits in Eq. Q and 
below: the limit uj — » is always taken in the dynamical 
sense, i.e. vq <C oj <C 2sf- Using Eq. Q in Eq. ([I]) and 
solving for uo we find that, to leading order in q, 

^pi(^^O) = Q • (3) 

In the same limit the imaginary-part of the low-frequency 
conductivity <j(uj) = ie 2 uJX P p(oo) /q 2 has the form 

f° 2 ?7 2 A 

9m a{w) -> . (4) 

UJ 

It then follows from a standard Kramers-Kronig analysis 
that the real-part of the conductivity has a ^-function 
peak at uj = 0: 5fte <j(uj) = V5(uj) where the Drude weight 

V = ne 2 v 2 A . (5) 

In the presence of disorder the ^-function peak is broad- 
ened into a Drude peak, but the Drude weight is pre- 
served. The Drude weight V defines an effective f-sum 
rule (cf. Ref. [21 ) in the dynamical regime vq <C uj <C 
2e F |221. 

We thus see from Eqs. (|3| and ([5| that the quantity A 
completely controls the plasmon dispersion at long wave- 
lengths and the Drude weight. In the following section we 
first relate A to the longitudinal pseudospin susceptibility 
and then carry out a self-consistent microscopic calcula- 
tion of A which demonstrates that its value is suppressed 
by electron-electron interactions. When this renormal- 
ization is neglected A ^ Ao and 

u; pl (q^0)=e F — -, (6) 



the RPA [23j (24J [25] result for the plasmon dispersion 
at long wavelengths. The RPA Drude weight V — 
^F^uni where, restoring Planck's constant for a moment, 
^uni = e 2 /(4h) is the so-called universal [26| [27^ l28| [29] 
frequency-independent interband conductivity of a neu- 
tral graphene sheet. We show that both uj p \(q — > 0) and 
V are substantially altered by electron-electron interac- 
tions. 



PSEUDOSPIN RESPONSE AND DIRAC-MODEL 
PLASMONS 

By using the equation of motion for the density 
operators [6 which appear in the density-density re- 
sponse function, Xpp(Q^) can be reexpressed in terms 
of the longitudinal current-current response function. 
When this procedure is applied to a Galilean- invariant 
system with mass m it leads immediately to the well- 
known result lim lim 5fte Xpp(q^) — nq 2 / (muj 2 ) . In the 

case of graphene, however, the current operator (defined 
as the derivative of the Hamiltonian with respect to fe) is 
directly proportional to the pseudospin operator [30] and 
we obtain instead 

2 2 

vq v q 

where & q is the component of the pseudospin fluctuation 
operator along the direction of qr, which we assume to 
be the x direction, and Xa x a x (q, w) is the longitudinal 
pseudospin-pseudospin response function. The latter de- 
scribes the response of a q to a pseudomagnetic field B q 
which enters the Hamiltonian with a term of the form 
G x _ q B q (notice that this has the opposite sign compared 
to the usual Zeeman coupling). 

Because of the presence of the infinite sea of nega- 
tive energy states Eq. ([7]) must be handled with great 
care [31] . In the noninteracting case we have 

([°q,P-q])= Yl [ COS (^) n i - - COs((p k - q )n^l q _] , 
\k\<Ak F 

(8) 

where cfk is the angle between k and the x axis and 
n k°^- = 1 is the occupation of the lower band. For 
q <C A/cp this can be rewritten as the sum of cos(<^k) 
over the region comprised between the circles |fe| < A/cp 
and \k — q\ < Ak^: only states deep in the negative 
energy Dirac sea contribute. When interactions are in- 
cluded is replaced by exact occupation numbers and 
additional terms associated with pseudospin-orientation 
fluctuations appear. In the next section we present a mi- 
croscopic time-dependent Hartree-Fock theory of A that 
is valid to first order in a ee . Since this theory neglects 
ground-state pseudospin and occupation-number fluctu- 
ations, which are of second order in a ee , the anomalous 
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commutator can be consistently evaluated in the nonin- 
teracting electron ground state and we find that [21] 



r vq A v 2 q 2 ge max 

hm -([a q ,p_ q ]) = 



j2 ^V 2 



(9) 



As we explain below, g£ m ax/(47r?; 2 ) is the negative of 
the pseudospin susceptibility of a noninteracting undoped 
graphene sheet, i.e. 

f^ = -Iimlim» B xS52- («,«), (10) 

where x^l x (Qi u ) * s the pseudospin-pseudospin response 
function of the noninteracting undoped system. Combin- 
ing the two terms on the r.h.s. of Eq. ([7]) we arrive at 
the following expression for A: 



A = lim lim 5fte 



(ii) 



Thus A has a very clear physical meaning: it is the pseu- 
dospin susceptibility of the interacting system regular- 
ized by subtracting the pseudospin susceptibility of the 
reference noninteracting undoped system. 

On quite general grounds it is possible to express the 
fully interacting value of A in terms of a small set of di- 
mensionless parameters by adapting to doped graphene 
sheets the original macroscopic phenomenological theory 
of Landau [6 , which is usually applied to normal Fermi 
liquids. In what follows, however, we will present a mi- 
croscopic theory of A, which we believe to be accurate 
at weak coupling and which enables us to draw quanti- 
tative conclusions on the impact of electron-electron in- 
teractions on the plasmon dispersion and Drude weight 
of doped graphene sheets. 



SELF-CONSISTENT HARTREE-FOCK 
MEAN-FIELD THEORY OF THE PSEUDOSPIN 
SUSCEPTIBILITY 

We now proceed to a quantitative microscopic calcula- 
tion of the pseudospin susceptibility A which goes be- 
yond the RPA. Specifically, we will take into account 
exactly the self-consistent exchange field which accom- 
panies pseudospin polarization. To this end we set up 
the time-dependent Hartree-Fock (HF) theory of the re- 
sponse of the system to a uniform pseudomagnetic field 
£? ext = B ext x oriented along the x direction. 

It is important to realize that the zero-frequency limit 
of the uniform pseudospin susceptibility is singular in the 
following sense. When £? ext is truly time-independent 
then its effect is simply to shift the occupied states in 
fe-space while reorienting the pseudospins. Changes due 
to fe-state repopulation and pseudospin reorientation at 
a given k cancel each other as required by gauge invari- 
ance, since a constant B ext can be eliminated by a gauge 



transformation [33]. At finite frequency, however, and no 
matter how small the frequency, the states are unable to 
repopulate (fe is a constant of the motion in the presence 
of the perturbing field) and pseudospin reorientation is 
the only effect that is left. This leads to a finite regu- 
larized pseudospin response in the zero- frequency limit. 
And this is clearly the limit of interest here [34] . 

The HF theory is usually described as an approxi- 
mate factorization of the two-body interaction Hamilto- 
nian into a product of simpler one-body terms. In the 
present case, the total HF Hamiltonian can be written in 
the following physically transparent form: 

Hhf = ^lal S *M k ) + • B ( k Mw , (l 2 ) 

where the HF fields Bo(k) and B(k) are defined by 
r d 2 k' 

Bo(k) = -J j^V k . k ,f + (k') , (13) 



and 



B(k) = B ext x 



vk-J 



d 2 k' 



with f±(k) — (n^\ ± n^ u; _)/2, where the nj^ are 
noninteracting band occupation factors, and n(fe) = 
B(k) /\B(k)\ is the unit vector in the direction of B(k). 
We have taken the limit q — > in Eq. (14) by considering 



(o) 



;V k - k *f-(k')n(k'),(U) 



a spatially homogeneous external pseudomagnetic field 
applied along the x direction. According to the previ- 
ous discussion the occupation factors in fe-space are not 
affected by the perturbation. The second term in this 
equation is the band pseudomagnetic field (see Wd), and 
the last term is the exchange field. The Hamiltonian TYhf 
has two bands with energies e^(k) = B (k) ± \B(k)\. 
In time-dependent HF theory electrons respond to the 
external field and to the induced change in the exchange 
field. 

In the absence of the external field n(k) = /c, so that 
there is no total pseudospin polarization, and |-B(fe)| — > 
|£? eq (fc)| depends only on k = \k\ [16]. Our aim here 
is to calculate the additional exchange field that arises 
from the polarization of the pseudospin when £? ext ^ 
0, since this determines the exchange correction to the 
pseudospin susceptibility. 

For q = the response is due to vertical interband 
transitions at wavevectors with \k\ > k^; transitions with 
\k\ < &f are Pauli blocked. Because only the transverse 
component SBt (oc z x fe) of the pseudospin field con- 
tributes to these matrix elements, it is evident that only 
SBt is relevant to the calculation of the susceptibility. 
We find that 

5B T (k) = -[5B Ti i(k) sinOfc)] z x k , (15) 
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FIG. 2: The exchange contribution to the equilib- 
rium pseudomagnetic field and the induced transverse 
pseudospin field. The data shown in this figure (as well as 
those reported in Fig. [3} have been obtained using a Thomas- 
Fermi screened potential V q = 27re 2 /[e(ga ee &F + <?)]• The data 
labeled by the blue solid line refer to the quantity Ea(x) plot- 
ted as a function of x = k/k? (for wavevectors k up to the 
ultraviolet cut-off /c ma x) for a ee = 0.2 and n^2x 10 13 cm -2 
(A = 20). The data labeled by the red solid line refer the 
solution u\(x) = 5Bt,i/ B ex t of Eq. (16) plotted as a function 
of x = k/kp for the same physical parameters. The dashed 
vertical line indicates the point k = /cf • Note that u\ is larger 
than unity, thus giving an induced transverse pseudospin field 
SBt,i that is larger than the bare external field B ex t- 



where 5B^^{k) / B ext = u A solves the integral equation: 
u A (x) = 1 + J dx'K(x,x f ) u A (x f ) (16) 

with 

Here we have introduced the dimensionless Coulomb 
pseudopotentials 



2tt 



dO 



2?re 2 J 2tt 



r <l\q=\k-k'\ ' 



being the angle between k and k' , and all wavevectors 
have been scaled with kp. The quantity 



\ r _ 

^a(x) = -a ee J dx'x'Vi(x, x') 



(19) 



is the exchange contribution to the equilibrium pseudo- 
magnetic field B eq (k). Illustrative numerical results for 
Ea and u A [the latter as obtained from the self-consistent 
solution of Eq. ( [l6| )] for a ee = 0.2 and A = 20 are shown 
in Fig. [2] 



FIG. 3: The uniform pseudospin susceptibility of a 
doped graphene sheet. The data labeled by filled sym- 
bols refer to the Hartree-Fock value of the ratio A/Ao, as 
calculated from Eq. ( 20 ) , as a function of electron density n 
(in units of 10 12 cm" z ) for various values of graphene's fine- 
structure constant a ee - The dashed horizontal line represents 
the prediction of the RPA [23] El [25], for which A/Ao = 1 
for every value of n and a ee - 



After straightforward manipulations we arrive at the 
following HF expression for the ratio between A and its 
noninteracting value Aq in terms of Sa(#) and i/a- 



A_ 



1 



dx 



u A (x) 



l + Z A (x)/x 



1 



(20) 



A plot of the ratio A/ Aq as a function of electron density 
for various values of a ee has been reported in Fig.[3j From 
this plot we clearly see that A/Aq is substantially lower 
than unity. According to Eqs. (|3| and ([5| this implies a 
strong reduction of the plasmon frequency and the Drude 
weight. 

In order to explain this result we show in Fig. 4 the 
response of the pseudospins to a pseudomagnetic field 
in the +x direction. Notice that the pseudospins in the 
lower band are tilted away from the pseudomagnetic field, 
while those in the upper band are tilted towards the pseu- 
domagnetic field. Because there are many more particles 
in the lower band than in the upper band we see that 
the total pseudospin response is negative. It is only after 
subtracting the noninteracting undoped response that we 
obtain the positive quantity A. It should be evident from 
this description that any many-body effect that enhances 
the total pseudospin response will reduce the value of 
A, while any many-body effect that suppresses the total 
pseudospin response will increase the value of A. 

In the present theory interactions affect the value of A 
in two competing ways, (i) The exchange field enhances 
SBt,i relative to £? ext (u A = SBt,i/ B ext > 1), thus en- 
hancing the total pseudospin response. From what we 
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4b) 




will be observable in experiments of the kind discussed 
in Refs. [261 [27J EH [29] . 

Finally let us comment on the broader implications 
of our results. Effects similar to those described in this 
article are also expected in graphene bilayers and other 
few-layer systems. The lack of Galilean invariance should 
also affect the cyclotron resonance frequency and oscil- 
lator strength when the 2D sheet of graphene is placed 
in a perpendicular magnetic field. Undoubtedly much 
interesting physics, potentially useful for applications in 
optoelectronics, has still to be learned from the study of 
graphene and other non- Galilean invariant systems. 
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FIG. 4: Response of the pseudospins to a pseudomag- 
netic field. Panel 4a): Pseudospins in the high-energy band 
tilt towards the pseudomagnetic field B ext applied along the 
+x direction. Panel 4b): Pseudospins in the low-energy band 
tilt away from the pseudomagnetic field. The reason for this 
unusual response is easy to understand. Pseudospins in the 
lower band are in their ground state: because of the anoma- 
lous sign of the pseudospin-pseudomagnetic coupling men- 
tioned in the main body of this article, they are anti-aligned 
with the intrinsic Dirac band pseudomagnetic field. When an 
additional external pseudomagnetic field is applied they will 
simply tilt away from it to minimize the energy. The occupied 
states in the higher band however, due to the Pauli principle, 
have pseudospins which are not in their ground state and 
which are aligned with the Dirac band pseudomagnetic field. 
When an additional external pseudomagnetic field is applied 
they thus respond in an unusual way tilting towards -B ex t- 



have said above it follows that this effect gives a negative 
contribution to A. (ii) The exchange contribution Sa to 
the equilibrium pseudomagnetic field increases [16] the 
conduction valence band splitting, thus suppressing the 
total pseudospin response. From the previous discussion 
it follows that this effect gives a positive contribution to 
A. Our calculations indicate that effect (i) dominates, 
resulting in a net reduction of the value of A. Physi- 
cally, the enhancement of the total pseudospin response 
(and thus the reduction of A) is a consequence of the 
gain in exchange energy that occurs when pseudospins 
in the same band tilt together towards a common direc- 
tion. We believe that the effect discussed in this article 
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